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Abstract. The relaxation rate and temperature-dependent switching field curve 
of a spherical magnetic nanocluster are calculated by including the effect of surface 
anisotropy via an effective anisotropy model. In particular, it is shown that surface 
anisotropy may change the thermally activated magnetization reversal by more 
than an order of magnitude, and that temperature-dependent switching field 
curves noticeably deviate from the Stoner-Wohlfarth astroid. With recent and 
future /^-SQUID measurements in mind, we indicate how comparison of our results 
with experimental data on isolated clusters may allow one to obtain valuable 
information on surface anisotropy. 



1. Introduction 

Understanding the dynamics and the mechanisms of the magnetization reversal in 
magnetic nanoclusters is essential for technological applications especially in the area 
of data storage. For such applications, rather fine nanoclusters are needed to increase 
the areal density, within ranges allowed by the superparamagnetic effect. On the 
other hand, in such small systems the magnetic state at the surface differs in many 
respects from that in the bulk. As such, the surface/interface effects cannot be ignored 
in investigations of the dynamics of magnetic nanoclusters. Now, surface effects are 
local effects and therefore such nanoclusters have to be regarded, in principle, as 
crystallites of many atomic magnetic moments. The problem with this picture of a 
nanocluster is that the study of the dynamics becomes a tremendous task because of 
the large number of degrees of freedom that have to be dealt with when computing, 
for example, the magnetization reversal time of the nanocluster. A way of overcoming 
this difficulty has been provided by the development of an effective model [T| which 
maps the many-spin cluster onto a macrospin representing the net moment of the 
cluster and whose energy contains mixed uniaxial and cubic anisotropies. The sign and 
magnitude of the effective coefficients of these anisotropy contributions depend on the 
size and shape of the cluster, crystal structure of the underlying material, and physical 
parameters such as the exchange coupling and local anisotropy constants. This model 
obviously provides a compromise between i) the Stoner-Wohlfarth (SW) macrospin 
approach [see the review article [2] and references therein] based on coherent rotation 
of all atomic spins and making use of the Neel-Brown theory for the calculation of 
relaxation rates, and ii) the many-spin approach [3, 4 J employing the Langer's theory 
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of the decay of metastable states [5]. In fact, these two theories provide asymptotic 
relaxation rates for each elementary process, i.e., the transition from a given state to a 
more stable state through a saddle point, and at this level they are fully identical and 
yield the same results. Finally, it is worth mentioning that this effective model also 
makes it possible to take into account more easily inter-particle dipolar interactions 
and therefore to provide us with a basis for interpreting experimental observations, 
e.g., through the zero-field-cooled/field-cooled magnetization, AC susceptibility, FMR, 
etc. 

In the present work, we compute the magnetization reversal time and demonstrate 
the effects of temperature and surface anisotropy (SA) on the switching field of 
a spherical nanocluster. We also indicate how, by comparing our results with 
experiments, it may be possible to gain information concerning SA in isolated magnetic 
nanoclusters. 



2. Energy and relaxation rate 

According to the effective model [I], the energy of a magnetic nanocluster may be 
written as 

E = - f i s U-m-K 2 Vm 2 z + ^K 4 V ^ m 4 a (1) 

a— x,y,z 

where m a , a — x,y, z, with |m| = 1, are the components of the verse of the nanocluster 
net magnetic moment and fj, s its magnitude, and V the nanocluster volume. H is the 
magnetic field applied at some angle ip with respect to the uniaxial anisotropy (UA) 
axis z. As mentioned in the Introduction, the parameters K2 and K& are effective 
anisotropy constants whose sign and magnitude depend on the intrinsic features 
of the nanocluster: lattice structure, shape and size, the (ferromagnetic) exchange 
coupling constant J, and on-site core and surface anisotropies. For example, for a 
spherical cluster with no UA in the core, the cubic anisotropy (CA) contribution in 
Eq. ([1]) is induced by Neel's SA and we have K4 = nAf K^/zJ, where N is the 
number of atoms in the nanocluster, z the coordination number, K s the on-site SA 
constant, and k (~ 0.53465 here) a surface integral that depends on the intrinsic 
features mentioned previously and also on the SA model. In the sequel, we use the 
dimcnsionless parameter £ = Ka/Kv- 

The calculation of the reversal time may be accomplished in two ways. The first 
is analytical and provides asymptotic expressions for the relaxation rate corresponding 
to one escape route, i.e., a passage from an energy minimum (1) to a minimum (2) 
through a saddle point (0). For this one may apply Langer's approach [5] which yields 
the relaxation rate in the form [7] 

where Z and Z\ are the partition functions computed in the neighborhood of the 
saddle point (0) and the minimum (1), respectively, and f2o is the damped frequency of 
oscillations about the saddle point (0). It is given by the absolute value of the unique 
negative eigenvalue of the transition matrix obtained from the linearized Landau- 
Lifshitz equation around the saddle point (0). 
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The application of Eq. ^ to a single magnetic moment and a standard steepest 
descents estimate of Z and Z\ leads to Brown's relaxation rate formula [5] 
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Here a is the dimensionless damping constant, (3 = 1/ksT, ks is Boltzmann's 
constant, T the absolute temperature, tn = /?Ms(l + a2 )/(2«7) is the free diffusion 
time. The coefficients are the Taylor expansion coefficients of the energy ([T|) near 
a stationary point (J). In general, the energy possesses several minima and saddle 
points. The calculation of the overall relaxation rate (the inverse relaxation time of 
the magnetization) must therefore involve a suitable combination of the elementary 
rates provided by Eq. ([2]) or ([3]). The escape rate rm^>(2) in Eqs. ©J® is valid only 
in the intermediate to high damping limit, i.e., for a > 1. For a< 1 the relaxation 
rate should be estimated using a universal turnover formula [llj 



T =A(aS)T 



(1)^(2), 



(4) 



where A(aS) is the depopulation factor defined as 

I / - x ln{l-cxp[-aS'(A 2 + l/4)]} 



A(aS) = exp 



A 2 + 1/4 



dX 



(5) 



and S is the dimensionless action variable. This formula is universal in the sense that 
it provides an accurate description of the relaxation rate for all values of the damping 
parameter a 





Figure 1. Energyscape for f = 1/2 (left) and 5 (right) and zero field. 

In the absence of applied field, it is easy to analyze the energyscape [see Fig. 
[T] and obtain analytical expressions for the relaxation rate using Eq. ([2]) or ([3]). 
In the case of dominating UA, i.e., for £ < 1, there two minima at 6^ = 0,n, 
with the energy Eq — —a(2 — C)/2, and 4 saddle points at 9^ — tt/2; ip^ — 
7r/4, 37r/4, 57r/4, 7ir/4 whose energy is £q ^ — leading to the energy barrier 



Magnetization reversal of a nanocluster 



4 



Ae = — £o'™' ) = c(4 — C)/4, where a — j3K 2 V . As the CA contribution increases 
there first appear saddle points at the equator (see Fig. [!}, then the minima become 
flat for Q — 1, and when £ is further increased the UA minima split and become local 
maxima surrounded by 4 deep minima. For large £ > these minima are located 
along the directions [±1,±1,±1]. For £ > 1, the saddles between the regions with 
m z > and m z < and the corresponding energies are the same as for £ < 1, 
while there are 4 minima for m z > and 4 minima for m z < that are given 
by cos6>( m > = ±y / (C + 2)/3C, <p {m) = tt/4, 3tt/4, 5tt/4, 7tt/4, with the corresponding 

energy £g = er(C 2 — 2£ — 2)/6C leading to the energy barrier Ae = fJ (C + 2) 2 /12C- We 
note that there are also saddle points connecting the minima within each of the regions 
m z > and m z < but they are irrelevant in the calculation of the magnetization 
reversal time. 

Hence, in the absence of applied field, we compute the relaxation rate for each 
escape route and multiply the result by the symmetry factor (= 8) that takes into 
account the number of escape routes. This leads to the following result for the inverse 
of the reversal time r = r _1 



7T\/2 



C_l1 e -a (C+2) 2 /(i2C) for C> i 
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for C < 1) 
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A C(C + 2 ) + (2-C). 
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For £ ~ 1, the parabolic approximation at the minimum breaks down, so that 
Eq. ((2]) must be used instead of Eq. ([3]) as Z\ can be computed beyond the parabolic 
approximation. Furthermore, it is clear that Eq. ^ cannot be used when ( — ► and 
an appropriate expression for this crossover to the pure UA has to be used [B]. 




Figure 2. Relaxation rate against £ for a = 20 in zero field. 

The behavior of the relaxation rate in zero field as a function of £ is plotted in 
Fig. [2] for a — 20 and a — 1. The full line is a plot of the formula in Eq. with Z\ 
being computed numerically in the whole upper plane of the energyscape [see Fig. [Tj. 
The triangles represent the result of the MCF calculation. 
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When a magnetic field is applied in an arbitrary direction, as is necessary for 
the calculation of the switching field, the loci of the stationary points must be found 
numerically. Accordingly, the relaxation rate can also be calculated numerically using 
the matrix-continued fraction (MCF) method as suggested by Risken [9] and developed 
in Ref. [10]. According to this method, the solution of the Fokker-Planck equation for 
the distribution function W($, cf>, t) of magnetization orientations may be turned into 
an eigenvalue problem for which the smallest non vanishing eigenvalue Ai yields the 
relaxation time r as r = A^ 1 . This relaxation time represents the time of reversal of 
the magnetization when the ratio £ is small, e.g., £ < 0.3. The MCF method provides 
a systematic check of the asymptotic formulae of the overall escape rate and allows 
us to bypass the energyscape analysis, which can become rather involved. On the 
other hand, the analytic calculation can be useful when the barrier is so high that 
the computation of the relaxation rate by the MCF method becomes extremely time 
consuming. 

As can be seen, good agreement is obtained for all values of (. One may see that 
the relaxation rate reaches a maximum that corresponds to a £ value £ max . This is 
due to the fact that for £ < £ max the barrier is lowered owing to the creation of saddle 
points, whereas for £ > £max the barrier is raised owing to a deepening of the minima. 
As a consequence the escape rate shows a bell-like shape. 

3. Surface and temperature effects on the switching field 

By definition, the switching field (SF) is the field Hg at which the magnetization 
overcomes the energy barrier and reverses its direction. At zero temperature, this 
reversal is only possible if the energy barrier is fully suppressed. This requires the 
application of an external magnetic field, called the critical field or the field at the 
limit of metastability. This definition corresponds to the SW critical field whose 
extremity describes the SW astroid, i.e., the limit-of-metastability curve [2]. At finite 
temperatures, the SF is intuitively smaller and relaxation of the magnetization must 
be accounted for. Moreover, experimental observation of the magnetization reversal 
depends on the relaxation time of the cluster and on the measuring time r m of 
the experimental setup. Therefore the magnetization reversal can be experimentally 
observed only if the relaxation time is in the time window of the experiment, or 
equivalently, if the relaxation rate is equal to the measuring frequency v m = t~ . 
Therefore, for the experimental observation of the magnetization reversal at finite 
temperatures, we must have 

r(<7, (,h s ,ijj,a) = v m (7) 

where hs = /i s Hg /(2K2V) is the switching field to anisotropy field ratio. This 
equation can be solved numerically for hs (c, £> ip> oi) and v m for given values of a, 
C, a and tp. At very low temperatures (a — * 00) and ( = this yields the SW astroid 
[5] ■ For the given parameters a, £, a, v m , hg can be plotted against the angle ip. For 
small C values, a parametric plot of the components h x = hs s'mip, h z = hgcosip, 
with hg obtained from Eq. |(7J) may be used. The field-dependent reversal time r is 
computed using the MCF technique for different values of the parameters a ~ T^ 1 
and C and the result is plotted in Fig. [3] In fact, we calculate the dimensionless 
relaxation rate Ttn with tn = t s x a, where t s = 1/ (o^^) with H2 = 2K2V/ 
determines the time scale within the system. For a spherical cobalt cluster of 3 nm in 
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Figure 3. 2D switching field for different values of a (~ T — 1 ) and f (cubic 
anisotropy contribution). The magnetic field is applied in the xz plane and the 
uniaxial anisotropy axis lies in the z direction. 

diameter with a = 1.0, t s evaluates to r s ~ 1.7636 x 10 s. In Fig. [3jthe measuring 
frequency has been taken as v m = 100 Hz. 

Two observations can be made from the results in Fig. [3] First, we observe a 
nearly homogeneous shrinking of the SF curve due to thermal fluctuations, and this 
has been observed in cobalt nanoclusters in Ref. [H] and also obtained numerically 
in Ref. [13j within the Langevin approach. The SF decreases with temperature owing 
to a decrease in the relaxation time. Moreover, Fig. [3] shows that upon increasing the 
CA contribution, i.e., when £ goes from to 0.25, there is a flattening of the SF curve, 
that is the critical field is reduced in the longitudinal direction (jjj = 0) and increased 
in the transversal direction (ip = tt/2). Recalling that £ (or K4) is related to the SA, 
these results show how the latter affects the switching field. This flattening of the 
SF curve owing to SA was observed in Ref. [14] where such a curve was computed at 
very low temperature for a spherical many-spin particle with local (on-site) anisotropy, 
uniaxial in the core and transverse on the surface with a variable constant. 

In the present work, we find that when the SA changes in intensity, the relaxation 
rate changes by more than an order of magnitude. More precisely, Fig. [1] shows that 
when £ = K4/K2 increases from (nearly) zero (dominating UA) and exceeds 2, the 
(dimensionless) relaxation rate Ttn increases from its minimum (~ 10 -6 ) and exceeds 
10~ 5 . Moreover, we note that since the relaxation rate has a maximum, optimum 
stability of the cluster magnetization against thermally-activated reversal is better 
achieved either at very small or large £. Indeed, for very small £ (dominating UA) 
there are no saddle points, and for large C (dominating CA) the energy barriers are 
very high, and thereby the relaxation rate, or the switching probability, is very small. 
In this context, these results imply that one may tune the thermal stability of the 
nanocluster magnetization by tuning its SA contribution. 
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These effects can be observed in /i-SQUID measurements on isolated (Co or CoPt) 
clusters especially at low temperature, and as a matter of fact, some preliminary 
measurements have indeed confirmed the flattening of the switching field curve but 
more measurements are needed to confirm this result. 

4. Conclusions 

The magnetization reversal time of a magnetic nanocluster is calculated taking 
account of surface anisotropy by using an effective macrospin model pQ that maps 
the energyscape of the many-spin nanocluster onto that of a macrospin with uniaxial 
and cubic anisotropics. The mapping allows us to compute the relaxation rate of a 
magnetic moment in the effective anisotropy potential as a function of the cubic-to- 
uniaxial anisotropy ratio so that surface anisotropy is accounted for. The analytical 
asymptotic formulae for the reversal time are given in both limits of weak and strong 
cubic anisotropy in the absence of applied field. The cubic anisotropy first creates 
saddle points, thus lowering the energy barrier, and then by further embedding of the 
minima it raises the energy barrier. The matrix-continued fractions method is used in 
order to check the range of validity of the asymptotic expressions of the relaxation rate 
in zero field. This method is also used to compute the reversal time for an arbitrary 
magnitude and direction of the field in order to calculate the switching field curves. 

These developments allow us to compute the switching field curves for weak 
cubic anisotropy as a function of temperature and intensity of the cubic anisotropy 
contribution. The switching field curve exhibits two main features. First, this curve 
homogeneously shrinks with increasing temperature and disappears completely at the 
blocking temperature of the nanocluster. This is in a qualitative agreement with the 
experimental observations on cobalt nanoclusters. Second, at a given temperature, 
we observe a flattening of the switching field curve when the intensity of the cubic 
anisotropy increases. More precisely, the longitudinal component of the switching 
field decreases whereas its transverse component increases. The more involved case of 
strong cubic anisotropy is left for a future investigation. Finally, we emphasize that a 
comparison between the present results and /i-SQUID measurements would provide us 
with an approach of further investigating surface anisotropy in magnetic nanoclusters. 

We thank William T. Coffey for reading the manuscript and making useful 
comments. HK acknowledges instructive discussions with Veronique Dupuis about 
the /i-SQUID measurements. 
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